Nonlinear polarisation and dissipative correspondence between 
low frequency fluid and gyrofluid equations 

Bruce D. ScottQ 
Max-Planck-Institut fur Plasmaphysik, 
Euratom Association, D-85748 Garching, Germany 
(Dated: February 2, 2008) 

Abstract 

The correspondence between gyrofluid and low frequency fluid equations is examined. The lowest 
order conservative effects in ExB advection, parallel dynamics, and curvature match trivially. The 
principal concerns are polarisation fluxes, and dissipative parallel viscosity and parallel heat fluxes. 
The emergence of the polarisation heat flux in the fluid model and its contribution to the energy 
theorem is reviewed. It is shown that gyroviscosity and the polarisation fluxes are matched by the 
finite gyroradius corrections to advection in the long wavelength limit, provided that the differences 
between gyrocenter and particle representations is taken into account. The dissipative parallel 
viscosity is matched by the residual thermal anisotropy in the gyrofluid model in the collision 
dominated limit. The dissipative parallel heat flux is matched by the gyrofluid parallel heat flux 
variables in the collision dominated limit. Hence, the gyrofluid equations are a complete superset 
of the low frequency fluid equations. 
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I. INTRODUCTION 



Low frequency reduced fluid equations are used to treat a variety of phenomena in mag- 
netised plasma dynamics, including turbulence HQ and tearing modes J3|. The usual 
derivation path is to start with the Braginskii collisional fluid equations 4|, and then to 
solve for the velocity in the Lorentz force rather than the inertia. To lowest order in the in- 
ertia/gyrofrequency ratio, balance among the principal forces is assumed, with pressure and 
the electric forces balancing the magnetic force. The result is a combination of E-cross-B 
and diamagnetic flow terms, arising from the electric field and pressure gradient, respec- 
tively. In the conventional magnetohydrodynamic limit (MHD: considering a single velocity 
for all species and a single, total pressure) the diamagnetic velocity is necessarily ordered 
small, but in general the electron parallel dynamics holds the electric and pressure forces 
to similar level. The latter is the "adiabatic response" which couples the electron pressure 
to the electric parallel current through pressure forces and compressional motion parallel to 
the background magnetic field. Hence the MHD ordering cannot be taken, and E-cross-B 
and diamagnetic flows are at similar level. 

The correction due to the inertia becomes the polarisation drift, which is so called because 
it is opposite for electrons and ions. The main contribution in MHD ordering is due to the 
time dependence of the electric field. Nonlinear advection of the velocity field represents 
the polarisation nonlinearity. It is responsible for maintaining drift wave self sustained 
turbulence 0, If], 0] and also any Reynolds stress flow phenomena si, S]. In general the 
diamagnetic contributions enter at the same order and give rise to what is generally called 
"gyroviscosity" — the cancellation of advection by the diamagnetic velocity in the equation 
of motion 10j. Once it is established that the polarisation drift enters at all, it is necessary 
to keep it also in the ion temperature dynamics, since the dynamics of the temperature and 
density are at similar order. The logical chain to this starts with the adiabatic response, 
including the compression in the electron density equation, then noting the equality of 
the electron and ion densities and the similar order in compression of the polarisation and 
parallel currents, and finally the polarisation drift entering the ion density and temperature 
at similar order 11]. 

Less familiar is the same phenomenon concerning the heat flux. The Braginskii model 
starts with a drifting Maxwellian distribution, with not only arbitrary velocity but also 
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velocity gradient, to lowest order in the inertia/collision frequency ratio. At next order, 
velocity gradients appear but heat flux gradients do not. This effectively and implicitly 
assumes that heat fluxes are subthermal: the heat flux is assumed to be smaller than the 
pressure times the velocity (see also Ref. [3] for similar considerations regarding implicit 
assumptions on the electron inertia and the magnetic current). However, this is not true even 
in the diamagnetic flows and heat fluxes, which in the presence of temperature gradients are 
of similar strength. The implicit assumption of small heat fluxes breaks down completely. 
When the MHD velocity ordering also breaks down, it follows that the diamagnetic heat flux 
is of similar magnitude as the pressure times the E-cross-B velocity. This has been noted 
before, by a treatment showing that the heat fluxes must be kept in the gyroviscosity even 
to obtain the standard fonn of the poh— eurrent Q. However , one has to go further 
and consider inertia in the formulation of the perpendicular heat flux itself 14j. This is one 



order higher in the moment hierarchy considered by Braginskii, which is why it is rarely 
considered. Nevertheless, polarisation enters the heat flux equation as the correction due to 
finite inertia upon the diamagnetic heat flux balance. Then, since the polarisation enters 
the density and temperature equations at the same order, and the polarisation heat flux 
and velocities are also of the same order, the polarisation heat flux should be considered in 
the temperature equation. We will review this herein as a preparation for establishing the 
correspondence between the gyrofluid and low frequency Braginskii equations. Ultimately, 
correspondence is found in the nonlinear advection effects only if the polarisation heat flux 
is kept in the fluid model. 

These polarisation phenomena enter the gyrofluid equations differently. The gyrofluid 
equations have an entirely different derivation path Jj], [3| , starting with t he g yrokinetic 



equation with the low frequency and small amplitude orderings already taken [17|. Polarisa- 
tion enters the charge balance equation rather than the density and temperature equations, 
since the latter are for the gyrocenters and not the particles themselves. The polarisation 
density balances differences in the gyrocenter densities, maintaining quasineutrality fl8 |. 
The time derivative of this gyrokinetic polarisation equation gives a relation analogous to 
the current balance (equivalently, vorticity) equation in the fluid models, with the time 
derivative of the polarisation density being the same as the divergence of the polarisation 
current. Underlying this is the Lie transformation between particle and gyrocenter coordi- 
nates at the gyrokinetic level Jjj] . Moments over this transform give the equations describing 
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the particle and gyrocenter representations of the moment variables, corresponding to the 
fluid and gyrofluid models, respectively. 

The gyrofluid equations are of significance because they allow treatment of this drift dy- 
namics at arbitrary order in the finite gyroradius parameter (generally, the square of the 
perpendicular wavenumber normalised to the gyroradius). Tearing modes and reconnection 
involve inertial layers which are thinner than the ion gyroradius 20|. Tokamak edge tur- 



bulence has a vorticity spectrum which always reaches down below the ion gyroradius [21]. 
Treatment of these is generally beyond the limits of equations whose derivation assumes the 
gyroradii are all small. Nevertheless, the low frequency Braginskii equations have a sys- 
tematic derivation, and it is desirable to know whether the gyrofluid equations correspond 
properly to these under the limits within which the Braginskii equations are perfectly valid. 
That task is the purpose of this work. In the linear MHD limit the correspondence between 
gyroviscosity and the finite gyroradius corrections in the ion density equation were already 
shown 151 ] . Examination of the nonlinear gyroviscous "force" in the MHD limit found cer- 
tain correspondences 22J. Herein, we complete the correspondence in the fully two-fluid 
limit. It is recovered only if the polarisation heat flux is kept in the fluid model. Viewed 
another way, this effect has always been present in the version of the gyrofluid model which 
keeps perpendicular and parallel temperature moments Q • The correspondence ques- 
tion is completed by examining the dissipation model in the gyrofluid equations concerning 
viscosity and parallel heat fluxes. Ultimately, the gyrofluid equations are found to recover 
the low frequency Braginskii equations, in the Braginskii limits of long wavelengths, small 
heat fluxes, and complete collisional dominance. 

The following sections respectively concern (II) the polarisation heat flux and its effect 
on the free energy theorem within the low frequency fluid equations, then (III) the cor- 
respondences concerning polarisation in the density and temperature equations including 
all the finite gyroradius nonlinearities, then (IV) the collisional viscosity effects including 
correspondence to the anisotropic corrections sometimes included in turbulence equations, 
and also the contribution of heat fluxes to the viscosity, and then (V) the parallel heat flux 
effects, whose correspondence is the easiest to show. The gyrofluid equations in question are 
from the most general GEM (Gyrofluid ElectroMagnetic) model 23|. They will be intro- 
duced piece by piece as needed. The collisional fluid equations are much better known - 
see, e.g, the recent model including the anisotropy effect in viscosity and the full polarisation 
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velocity treatment in Ref. 



241 ]. A concluding commentary section (VI) is given at the end. 



II. POLARISATION INCLUDING THE HEAT FLUX 



Low frequency fluid equations can be derived directly using the equation of motion for 
each species 
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solving for u in the Lorentz force term to lowest order in oj/Vt c <C 1 and assuming an 
electrostatic perpendicular electric field with potential (justified by u <C k±VA] i.e., the 
dynamics is too slow for dynamical Alfvenic compression), 
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noting this gives solely the perpendicular component. The parallel component has its own 
equation, derived separately. The polarisation corrections are found by inserting this ui^ 
form into the inertia terms, 
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assuming flute mode ordering wherein WiiVii < ui ■ V. This is the standard version 10], 
usually behind the derivation of the equations in turbulence models. Alternatively, a system- 



atic procedure splitting the velocity into solenoidal and parallel pieces [25(, which can also 
include a potential- flow compressional piece [26j, may be used. The solenoidal flow potential 
jecomes under MHD ordering or generally a combination of <fi and p. If drift ordering 
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is then taken, the equations become identical to the reduced forms. Drift order- 
ing refers to the small amplitude but unity-order nonlinearity limit used in the turbulence 
models. 

Application of drift ordering to the velocity and including the diamagnetic pieces in II 



results in the following form 
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where the d/dt operator includes the nonlinear E-cross-B advection 
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The last term in Eq. (jlj) is the polarisation velocity, whose charge-flux divergence is given 
by 

V., iZe u p = -V.-^-(v^ + — V XP ) (6) 

In the conventional gyro-Bohm normalisation for a single component plasma with singly 
charged ions this becomes 

V • u p = -V • ^V ± (0 + Pi) = -V • j f V±W (7) 

where Pi = r^nj + Tj) and Tj is the background ion/electron temperature ratio. The pressure 
gradient is linearised, and each species has its density and temperature normalised to its 
own background. The flux and velocity divergence enter the same way because of the 
normalisation and ordering. The total ion force potential is denoted as W. The time scale 
inferred by the divergence of the velocity is normalised to the profile scale and the sound 
speed, L±/c s , where c 2 = T e /Mi and L Te is usually used for L±. The double perpendicular 
derivative is normalised to the square of the drift scale, p 2 = c 2 T e Mi/e 2 B 2 , and it is useful 
to note that r^p 2 = p 2 , which makes the role of the finite gyroradius explicit. The above 
considerations constitute what is also called local ordering in the context of turbulence 
computation. 

A similar treatment for the heat flux starts with the (Mt> 2 /2)v moment of the kinetic 
equation, analogous to the Mv moment and the equation of motion. The heat flux equation 
is given by 

^ + ^VT = ^qxB (8) 

dtZ,2M Mc V ; 



under drift ordering (Eq. 11 of Ref. [I4j, after the diamagnetic cancellation is taken), Solving 
this to lowest order neglecting the inertial effects yields the diamagnetic heat flux, 

Using the same ordering and normalisation conventions as for the velocity, we find the 
divergence of the polarisation correction for singly charged ions, 

V ■ q p = -^r,V ■ J^±Ti (10) 

under gyro-Bohm normalisation (q is normalised the same way as pu). Clearly, this is the 
same order as the velocity polarisation divergence if the gradients of the state variables (po- 
tential, densities and temperatures) are all comparable. The polarisation velocity divergence 



enters both the density and temperature equations, and the polarisation heat flux divergence 
enters the temperature equation. This sort of consistency is well known for the diamagnetic 
fluxes themselves, since with similar (e/T)V0 and VlogT and Vlogn they all enter at the 
same order with each other and with toroidal compression of the E-cross-B velocity {29! . [3^] 
(cf. discussion and manipulations in Ref. 1 3 11). but the same results concerning the polari- 
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sation fluxes was not widely known before Refs. [131 . 1141 ] , and is still routinely missed by low 
frequency fluid models. 

Under the above considerations we have the normalised ion density and temperature 
equations with polarisation divergences, 

^1 _ V ■ ^VW + V||«|| = K{W + G) (11) 

3 rIT- d ^ d ^ 

2li ~ V ' Jt VW ~ 2 TlV ' d~t VTl + V||(?i|1 + qill) = K{W + G) + 2 nm) (12) 

The parallel velocity divergence is included in both equations and the heat flux divergence 
in the temperature equation. The terms denoted by /C are the remnant divergences of the E- 
cross-B and diamagnetic velocities (represented in total by W) and diamagnetic specific heat 
flux (represented by 5rjTj/2), due to the inhomogeneous magnetic field, after the diamagnetic 



32] . The curvature operator is then 



cancellation is taken in the temperature equation 
defined, e.g., 

/C(0) = -V-v £ = -V--|jBxV0 (13) 
in terms of the E-cross-B divergence. The quantity G in the ion density and temperature 



equations arises from thermal anisotropy. It is given by 241 ]. 



G 



0.96 
12^ 



K(W) - 4Vpn (14) 



in the collisional limit and represents viscous dissipation, with Vi the ion collision frequency 
normalised to c s /L±. 

For the electrons the convention is to neglect the mass everywhere except in parallel 
inertia (entering the parallel velocity and heat flux equations). The electron density equation 
is given by 

^ + V||«n=/C(0-p.) (15) 

in which anisotropy and polarisation (electron viscosity and inertia) are neglected. The 
electron pressure gradient is linearised in the same was as for the ions, with the minus sign 



reflecting the normalised temperature/charge ratio. The quasineutrality condition is given 
by the subtraction of the two density equations and neglecting the space charge density, so 
that 

V • jVW = V||(«n - u||) - K{ Ve +Pi + G) (16) 

equivalently, V • J = 0, whose three pieces are the polarisation, parallel, and diamagnetic 
divergences, respectively. We note that 

J,, = u|| - u|| (17) 

defines the parallel current (under the normalisation); this is usually used to eliminate v\\ in 
favour of J\\. We may further subtract this from the ion temperature equation to obtain 

Q A r P. K A K 

2 ~d ~ 2 nV ■ j t VTi + V|1 + %||) = w - Pe) + r MTi) (18) 

eliminating the polarisation divergence in the velocity but not the heat flux. In this equation, 
the explicit ion velocity divergences are replaced by the electron ones, but the ion heat flux 
divergences remain. These are the polarisation and diamagnetic heat flux terms, respectively 
the second and last terms in the line above. 



A. Free energy in the fluid model 

The complete set of equations in the fluid model is given by 

V • jVW = V,| J,, - /C(p e +Pi + G) (19) 

ATI 

^ + V| R |=/C(0-p e ) (20) 

YdT + V||(w " + ge||) = W " Pe) ~ \ K{Te) (21) 

3 dT- 5 d 5 

2^t " 2 TlW ' Jt VTi + V||(U|1 + = K ^ ~ Pe) + 2 nKiTl) (22) 
diii II 

^ + V||(p e + Pl + 4GO = (23) 

+ + V||(0 -p e ) = -i? el (24) 

where and n e are equivalent, v\\ is given by u\\ — J\\, and the parallel heat fluxes q e \\ and qq 
and the resistive dissipation R ei are left undetermined (at this level they may be given their 



Braginskii dissipative formulae 4j). The factor of 4G in the parallel momentum equation 
is also the result of anisotropy. Except for the retentions of the polarisation heat flux in 
Eq. 
Ref. 



and the electron inertia in Eq. (I24p . these equations are the same as those given in 
24J . The normalisation convention is the standard gyro-Bohm one, with Vn normalised 
against L±, not qR, which is why the un-scaled forms for 

Airp e 



0e 



111,, 



(25) 
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are used. The factor n gives the background Tj/T e ratio. The pressures are linearised as 
above. 

The free energy of the system is given by 
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where / dV denotes complete spatial integration and now the 7$ factors are put in explicitly, 
xcept for the last term, due to the polarisation heat flux, this has been analysed before 



33l . l34j . Insertion of Eqs. (I19ti24p into dS/dt finds this time derivative to vanish except for 



the dissipative terms (here, there are no gradient source terms since the profile gradients are 
kept within the dependent variables of the model; without explicit sources this corresponds 
to decaying cases initialised with a finite profile and a random bath of fluctuations, as in 
Ref. [35] ) . 

The last term in Eq. (1261) represents the polarisation heat flux. If Eq. (1221) is multiplied by 
TjTj and integrated, the d/dt terms yield the two terms in Eq. ([26]) explicitly dependent upon 
Tj. If Eq. ( IT9]) is multiplied by W and integrated, the d/dt term yields the term in Eq. (126]) 
explicitly dependent upon W. The resulting term W"V\\ J\\ is balanced by the contributions 
Tjn e and TjTj times Vy J\\ in Eqs. ( l20][22l and the contribution JmV||0 coming from Eq. (T24l) . 
And henceforth. The appearance of a heat flux term in the energy may be unfamiliar in a 
fluid model, but it is known from the gyrofluid model (cf. Ref. 36j and below) and also from 
treatments of extended fluid dynamics 37]]. Due to the relation between the polarisation 
heat flux (as, up to coefficients, the time derivative of the curl of) the diamagnetic heat flux, 
the energy contribution is equivalent to 2/5 times the square of the diamagnetic heat flux 
(up to normalisation). This is equivalent to the relation between the polarisation velocity 
and the E-cross-B and diamagnetic velocities, and the appearance of the square of the latter 



in the form of the perpendicular kinetic energy. For general gradient driven turbulence, the 
heat flux and velocity pieces are of comparable magnitude. 

The rest of the fluid model is given by the dissipation due to parallel heat fluxes and 
viscosity. The latter comes from temperature anisotropy, in this case G. Starting with a 
diagonal pressure tensor 

P = diag{p ± ,p ± ,p||} (27) 
we split it into isotropic and traceless parts, 

P = pg + n n = 2Gdiag{-l,-l,2} (28) 

where g is the metric tensor, identifying 6G with the anisotropy Ap = p« —p± (see Section IV, 
below). With the factor of density common, Ap is equivalent to AT. The viscosity model sets 
collisional dissipation of G against perpendicular and parallel velocity divergences as above. 
Similarly, the parallel heat flux formulae are given by setting their collisional dissipation 
against the corresponding parallel temperature gradients, e.g., 

r)/2 mq = —Ti^Ti (29) 



t> it>\\ rj 



4| with Ki = 3.9 results. 



for ions, with the coefficient set such that the familiar formula 

The remainder of this paper is concerned with recovery of these formulae (polarisation 
effects in the vorticity and ion temperature equations, viscosity through anisotropy, and the 
parallel heat fluxes), from the gyrofluid model under the same ordering conventions as for 
this one. 

B. Interlude — bracket notation 

In several treatments of the equations of turbulence in confined plasmas the nonlinearities 
are explicitly written in a form which makes their conservation properties obvious. Basically, 
e.g., we ■ Vn e is written as [</>,n e ], where the bracket involves the perpendicular derivatives 
involved in the drift motion. It is variously written as 

dx dy dx dy 

in slab or local fluxtube treatments, or as 

I ( <90<9n e dn e d4>\ 

^ ne] = r{d^W-^de) (31) 
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in "cylinder" treatments, which with the field aligning coordinate transformations x = r 2 /a 2 
and i)k = q{6 — Ok) — C and s = 9 with q = q(r) becomes 



[0, n e ] = -7 



<9x \ dy 



dn P dn 



dn 



ds 



e I ^0 <90 



(32) 



at s = 6^ since dy^jdr vanishes there. Under fluxtube ordering d/ds is small compared to 
either d/dx or d/dy and the factor of q is replaced by a constant, and Eq. f)32p reverts to 
the form in Eq. (|30|) with a multiplier of 2q/a 2 which can be normalised away. More detail 
on this is given in Ref. 3s|. 

With the bracket notation we make use of the following properties 



JdV[f, g ] = o f[f,g] = \[f 2 ,g] 



g[f,g] = ^if,g 2 } 



(33) 



in the manipulations below. That is, the bracket is a perfect divergence, and both energy 
and entropy are conserved. Useful manipulations include 



V1[/,0] = V-[/,V.l<7]+V.[V_l/,0] 
V ■[f,V ± g} = [V ± f,V ± g} + [f,Vlg] 



(34) 
(35) 



In any of these the perpendicular subscript may be regarded as understood, as is the con- 
traction implied by [V/, Vp]. 



C. Free energy, adiabatic response, and MHD ordering 

The fluid model's vorticity equation in this notation is 

j t V\W + [V ± 0, V ± W] = V,| J,, - K(p e +Pi + G) (36) 

where use is made of 

d d r , . ,„„. 

7t = M + ^ ] (37) 

Since the bracket is antisymmetric, the second term in Eq. (136]) is equivalent to [Vj_0, Vjj>j], 
and hence this "gyroviscous correction" is a proper warm-ion effect. The polarisation terms 
may be manipulated to show 

^y\w + v • [0, v ± w] = ^v 2 ± + v • [w, v ± <p} + vi^ (38) 
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Under MHD ordering the last term on the right hand side is dropped, as was done in Ref . 10] , 
unfortunately without the explicit statement that this depends strictly on V±pi <C n e eVj_</> 
remaining valid (un- normalised units). 

MHD ordering implies the lack of a rjVy Jy term in the energetics. This can only be 
reconciled if Vy Jy is neglected in all of the continuity equations, not just the one for n e , 
since quasineutrality ties the ion and electron dynamics together. If VyJy is neglected 
then V|ip e itself must be neglected in Eq. (1241) . Hence, we would be back not only to MHD 
ordering but to reduced MHD itself. This may be demonstrated by alternatively multiplying 
Eq. ( |T9l) by or W and integrating, and observing the logical consequences. Either reduced 
MHD is taken in its entirety, or the consequences of the adiabatic response are taken to 



their conclusion. No intermediate version is energetically closed 

It is the adiabatic response in the electron dynamics which disallows the use of MHD or- 
dering in this kind of turbulence. In Eq. (1241) the two static force terms (those not dependent 
on Jy or A\\) are often the largest, and even for edge turbulence there is a partial cancellation 
between them. Their difference determines Jy, mediated by induction, inertia, or resistivity 
according to whether finite /3 e , finite /i e , or finite v e is the strongest. Once Vyp e is kept 
in Eq. ff24|) then VyJy must be kept in Eq. (120]) . since in the energetics n e V\\J\\ + J\\V\\n e 
must become a total divergence. Since it is p e which appears in Eq. (124"]) . then VyJy must 
also be kept in Eq. ff2TT) . Now, the densities n e and are equivalent, so the appearance 
of VyJy in Eq. (|20l) implies the appearance of V ■ u p in the equation for rii, which is made 
explicit if Eqs. ( l20lfT9l) are subtracted (eliminating the VyJy term). To conserve against the 
implied TiVy Jy term, the pressure effects in polarisation must be kept since Eq. ( JT9l) must 
be multiplied by W rather than 0, since W contains 7in e . However, W also contains r^Tj. 
Hence, Vu p must also be kept in Eq. fl22|) . wherein it has been replaced by VyJy. Only now 
is the energetic loop started by Vyp e in Eq. (1241) closed, since now — WVy Jy closes against 
Ti(ni + Ti)V\\ Jy from Eqs. ( 12011221 as well as against — JyVy0 from Eq. ( !24l) . Hence the logical 
chain: adiabatic response in the equation for Jy, parallel compression in the equations for 
n e and T e , quasineutrality, by which VyJy implies V • u p , and finally the qualitative simi- 
larity among densities and temperatures. All this is forced by the similarity in magnitude 
among Vj_{0, n e , T e , Tj}, caused by the adiabatic response. Similar consequences may be 
found among the curvature terms — essentially, the retention of diamagnetic compression in 
the continuity equations also forces the complete two fluid version of ion polarisation. The 
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above analysis has been given for the local form of the equations; the corresponding one for 
the global form was given in Ref. 11]. 



D. Interlude — relation to extended fluid and extended MHD models 



The low frequency fluid models are themselves distinct 



dynamics or extended MHD 



39 



40 



41 



42 



43|, 



44, 



45 



rom what is called extended fluid 



461 ] . In the latter the low frequency 



is invoked to obtain expressions for heat fluxes and viscosities. These are not, however, 
expressed in terms of drifts and polarisation, but left in the native form with velocity and 
heat flux vectors. Specifically, the steps in Eqs. (jTH^J) are not taken. Moreover, explicit time 
dependence of the heat flux and dissipative viscosities are neglected (thermal anisotropy is 
assumed to be small - see below). Fluid drift theory replaces the vector forms with scalar 
quantities in the list of dependent variables; for example, and p and_U|| are the variables 
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47|], to be 



with which u is described. Further to that are the gyrofluid models 
discussed below. No matter the complexity, extended fluid models and fluid drift models 
Dreak down when k±pi becomes unity or larger, as it always does in tokamak edge turbulence 
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47 



48|, 



49j. Gyrofluid models are required to overcome this. Of course, since they have 



a different formulation, it is desired to know how well they recover the fluid forms when the 
latter are valid. That is the point of this work. 



III. GYROFLUID FLR NONLINEARITIES AND FLUID GYROVISCOSITY 

The gyrofluid model has a different structure from the fluid one — like the underlying 
gyrokinetic model, the moment variables (the model for the kinetic distribution function) 
are advanced independently for each species, and then the field equations (polarisation and 
induction) are solved for the electrostatic and parallel magnetic potentials. Analysis of the 
gyrofluid moment equations in the various limits proceeds the same way for each species. 
Due to the correspondences involved we concentrate mainly on the ions. Consideration of 
thermal forces at the end will then involve the electrons. This section is concerned with 
the gyroviscosity effects in the fluid model, meaning essentially all the differences in the 
polarisation between the general one and the MHD one (the latter involving only). We 
will show how these emerge naturally from the finite gyroradius (FLR) nonlinearities in the 
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ion gyrocenter density and temperature equations in the limit of small k±pi <C 1. 

The equations under consideration are for the density and the parallel and perpendicular 
temperatures (Eqs. 99,101,102 of Ref. [231), 



+ [<I>g, ni} + [Sla, T i± ] + V||«n =KU g + 



Pi\\ +Pi± + 



1 dT - 1 

2~m +2 [0G ' T4||] + V|l(M|l + gi|lll) = /c 



2v 4 G 



(39) 
(40) 



dJh 
dt 



+ [0G,T 4 . 



p G) n i + 2T iA ] + V\\q. 



tJ-l 



1 



/C 



6 G + 4fi G + p i± + 37*7} 



+ 2uiG 



(41) 



2 V 2 

where in terms of the perpendicular and parallel pressures the isotropic one is Pi = (p^ii + 
2]?i±)/3 and the difference is G = (py —p±)/6, dissipated in the term proportional to V{G. 
Normalisation is to a common background temperature T , with n = Ti/ZT giving the 
temperature/charge ratio. The factor of r, is folded into the pressures. The heat fluxes are 
also broken up into parallel transport of perpendicular and parallel energy, the (Mw^/2)w\\ 
and (Mwj_/2)wu moments, respectively, where w is the kinetic velocity in the co-moving 
reference frame. In the gyrokinetic and gyrofluid models is not used directly, but as w\, 
specifically, the magnetic moment \i = Mw±/2B, due to the low frequency ordering. These 
heat flux pieces are <fc|ii| and qi±\\, respectively. Finally, the potentials 0g and fie represent 
the FLR treatment. The Pade approximants are 

-b 2 <f)/2 



} G 



n 



1 + 6/2 " G (l + fe/2) 2 

in wavenumber space, with argument b = k\pf. In the limit of k\ — > 

- b -> rrfl <j) G -> (1 - 6/2)0 n G -> i 

which we will use to show correspondence. The notation of Ref. 



FLR treatment follows Refs. 



conservation as discussed in Ref. 




(42) 
we have 

-6/2)0 (43) 
2311 is used, and the 



with the necessary modifications to restore free energy 



A. The isothermal version 



To make the analysis easier to follow, we start with the isothermal model which neglects 
all considerations of temperature dynamics, including the anisotropy and heat fluxes. In 
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this case the ion and electron density equations are given by 

dm 



Of 

and they are related through the polarisation equation, 



+ [<f> G , rii] + V||U|| = K{(j> G + nm) (44) 
P + [0, n e ] + V||V|| =K,{(j>-n e ) (45) 



+ 7J- 1 — -^n. (46) 



1 + 6/2 4 1 + 6' 



using the Pade approximants (cf. Refs. 16j, |23j). This determines the gyrocenter density n 



in terms of the particle density (equal to n e ) and the polarisation contribution (due to 



as 



Ui = (1 + 6/2)n e + T r l b(j) (47) 

expanding in powers of 6 and keeping the 0(1) and 0(6) terms. This is now converted into 
configuration space identifying b with — TjVj_, leaving 



= n e - Vl0 - ^Vlpi (48) 



where Pi = Tin e . The FLR-corrected potential is given as 

<Pg= (1- 6/2)0 ^0 + |vi0 (49) 

up to 0(6). Here and below, is normalised against p~ 2 , so that pfV\ becomes TjV^. 

Now we use the equations for n e and rij to find the vorticity equation, using these forms 
to eliminate Ui and 0g i n terms of n e and <fi. First, the ion density equation becomes 

^ (n e - Vi0 - \^\p) + [0, n e ] - [0, Vi0] - i [0, V\ Pl ] + ± [V^, ft] 

+ V||U|| =JC{<f>+ Pi ) (50) 

where under K, only the 0(1) terms are kept. The terms under d/dt are from n.j. The first 
three bracket terms are from [0, rij], and the last bracket term is from the difference 4>g — <P- 
We manipulate the bracket terms involving as follows 

[4>, V\4>] = V • [4>, V ± 4>] - [V ± 0, V ± 0] (51) 
[(f), Vift] = V • [0, VxPi] - [V±0, Vxft] (52) 



15 



[V10, pi] = V ■ [V ± 0, ft] - [V±0, V±pi] = VI [0, Pi] - V ■ [0, Vxft] - [V ± 0, V ± p,] (53) 
so that with cancellations (noting also [Vj_0, V^0] vanishes) we obtain 

^ (n e - Vl0 - \^lP^j + [&n e \ - V • [0, V ± 0] - V • [0, V^] + ~Vi[0,ft] 

+ V||U|| = /C(0 + Pi ) (54) 

Then, under we replace the bracket [4>,pi\ with —dpi/dt noting all the other terms are 
0(b) corrections to the Vii and /C terms, so that 
d 

(n e - Vi0 - Vlpi) + [0, n e ] - V ■ [0, V ± 0] - V ■ [0, Vxft] 

+ V||U|| = /C(0 + pi) (55) 



<9f 



Finally, we combine d/dt and [0,] into d/dt, and and pi into obtaining 

^ - V • U/JV + V F || = /C(0 + ft ) (56) 
at at 

This is the ion density equation in the isothermal fluid model, wherein pi = T(n e and by 
quasineutrality the particle (not gyrocenter) densities are equal. Subtraction of Eq. (156"]) 
from Eq. (I2U|) in Section [II A[ we find 

V • jVW = V,| J,, - (1 + 7i)/C(n e ) (57) 

where we have inserted u\\ — fn = J\\ and p e +ft = (1 + Tj)n e . This is the same as Eq. ( |19p in 
Section III A\ under the isothermal gyro-Bohm normalised forms p e = n e and Pi = Tin e and 
G = 0. 

The rest of the isothermal fluid equations are 

dn 

-^ + V||(u||- J|,) = /C(0-n e ) (58) 

^ + (1 + T;) Vp e = (59) 

„ cM.ii dJii „ . , . . 

&~5r + ^ e ~u + V ll W ) ~ n ^ = -0.51/x e i/ e J|| 60 
at at 

and these satisfy the energetics as given in Section III At without the temperature dynamics 
(the latter inclusive of G and the heat fluxes). The resistive dissipation includes the 0.51 
coefficient from Ref. [J] and the collision frequency v e is normalised against c s /L±. What 
can be termed "gyroviscous correspondence" to the fluid model is thereby proved for the 
isothermal case, for any occurrence of nonlinearity within the local ordering. 
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B. With temperature dynamics 



Now we return to the version of the model with all the temperature dynamics, involving 
Tj|| and T^. We will make use of the particle representations of these variables. For the 
density the relation between the particle and gyrocenter representations is given by the 
polarisation equation, Eq. (I4U|) . whose thermal version is (Eq. 92 of Ref. [231] ). 

r -i 



Tin; + r 2 Tj. 



n,. 



Ti 



(61) 



The Ti and T 2 are the gyroaveraging operators. Eq. (161 p is the closure approximation to 
the gyrokinetic polarisation equation, 

J F M J -F M 



E dW 



eJ (8f) + e 2 



T 







(62) 



where 5f denotes the distribution function, / dW represents integration over velocity space, 
the sum is over species, Jo with argument k±v j_/fl acting upon 5f or is the orbit averaging 
operator, F M is the background Maxwellian, and e and T with n and M are the species con- 
stants giving the charge, background density and temperature, and mass. The polarisation 
equation comes from setting the particle charge density to zero 18], and the polarisation 
;erm itself ultimately arises from the transformation from gyrocenter to particle phase space 
ljj. The closure approximation to J dW F M J is Ti, with arg ument b. The closure form 
is given by T^ 2 by correspondence to linear kinetic theory [15(. The second operator r 2 is 
given by the logarithmic derivative of Ti with respect to b no matter the form chosen for 
r 1; since 



T 



or 



The Pade approximant forms for Ti and T 2 are 

1 



1 + 6/2 



1 ) F M J 



-b 2 /2 



(63) 



(64) 



as given in Ref 
Eq. (I42p result from 



16] . Hence, the gyroaveraged potential and its FLR correction as given in 



n 



G 



(65) 



That these are the same operators as the ones in Eq. (161j) . and that <pG and Qq are associated 
with rii and T^±, are fundamentals underlying the free energy conservation of the model 23]. 
The form of the nonlinear terms in the Tj^-equation results from the next higher moment 
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with jj,B and applying free energy conservation as a constraint. Further detail on this and 



FLR closure in general is given in Refs. 



23 



50|], which update Refs. [l5|,lR|. 



In these terms the particle (space) representations for the three state variables for the 
ions are 

r — i 

n sp = TtUi + r 2 Ti_L + — (66) 

Ti 

Ti\\ sp = TiTi\\ (67) 

T i±sp = T 1 T i± + T 2 ( ni + 2T i± ) + 2^0 (68) 

all arising from corresponding moments of Eq. fl62l) . The first is the same as the polarisation 
equation. In the second, resulting from the Mw^ — T moment, the parallel and perpendicular 
velocity space integrals separate, and the piece vanishes. In the third, the moment of 
/j,B — T over the Jq operator gives rise to the same factor of (I\ + 21^) as in the nonlinearities 
in the Tjii equation (Eq. |4~T1) itself, in addition to the term I^I^ coming from the moment of 
fiB — T over the Jq operator. 

As in the isothermal case, the fluid equations are found by constructing the time deriva- 
tives of these variables in the particle (not gyrocenter) representation. The \ow-k± limit is 
taken, with the 0(b) corrections kept only in the nonlinear advection terms. The equations 
for rii and Ti±_ occur together, ultimately due to the way J through its 6-dependence mixes 
the perpendicular moments. Up to 0(b) the particle representations are given by 

n isp = n i + V 2 ± (l)+^Vl(n i + T l± ) (69) 

T l±sp = T l± + Vi0 + |vi(ni + 3T l± ) (70) 

+ (71) 

with the terms on the right sides understood to be in the gyrocenter representation. In the 
0(b) terms the representations are equivalently the particle or gyrocenter ones, so that the 
inverses of Eqs. (16914711) are given by 



m- Vi0-£vi(ni + 2U) (72) 



2 



T i±gy = T i± - Vl0 - |vi(n, + 3T l± ) (73) 



Ti\\„ = ^|| - T i ViT.ii (74) 
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with the terms on the right sides understood to be in the particle representation. The 
gyroreduced potentials in Eq. fl65|) by 



; (T . = o + ^Vl0 fi G = |vi0 (75) 



The partial time derivatives of Eqs. fl72H74|) are taken, and then Eqs. f )39l - H0l) are used to 
evaluate the right hand sides, as was done in Eq. ( J501) above. 
For the density the result of the substitution is 

JL ( n . _ v 2 ± - -V^) + [<f>,n e ] - [<p, Vl<P] - ^[(f>,V 2 ±Pi ] + \{Vl^ Pl ] 

where V 2 ± Pi = TjV^(nj + Ti±), linearised as before, represents the combining of the V^rij 
and V\Ti± terms. The manipulations of the Vj_ operators are done exactly as before, and 
the result is 

£-V.!™ + V,., = K (, + ^) (77) 
For the perpendicular temperature the result of the substitution is 

j t (t i± - Vi0 - l - V 2 ±Pl - T& 2 ± T i± 

+ [<f>,T iX ] - [<f>, Vi0] - ^[<f>, V\vr] - [0, Vir^J + ^[Vl4>, Pl ] + [Vl0,7i2U] 

+ V||gijL|| = fC I \+2v t G (78) 

where FLR corrections to the /C terms are dropped as before. The nonlinear terms propor- 
tional to Ti± arise from the extra factors of 2Tjj_ in Eqs. ( HTlr73l) . The manipulations of the 
V;}_ operators are done exactly as before, and the result is 

~dt ' dt ~ ' dt V II*-MI = ^ I 2 j ^ ' 

with the last of the nonlinear time derivative terms representing the perpendicular part of 
the polarisation heat flux. For the parallel temperature the result of the substitution is 

| (r q - £v5r iB ) + [4>,T q ] - T p,vlT q ] + f [v^r,,,] 

+ 2V||(«|| + g i ||||)=/c(0 + p i |, + 2r i Ti|,)-4i/ < G (80) 
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where FLR corrections to the K, terms are dropped as before. The nonlinear terms propor- 
tional to Tj|| arise from the extra factors of TJn in Eqs. ( I40"f74l . The manipulations of the 
V\ operators are done exactly as before, and the result is 

2^r " 2 v ■ Jt v ^ + V|l(n " + ^ = K [ 2 J - 2 ^ G (81) 

with the last of the nonlinear time derivative terms representing the parallel part of the 
polarisation heat flux. The two temperature equations (Eqs. I79|8TI) are added to provide 
the final temperature equation 

If - v ■ - h v ■ s v - T > + v «<"« + *■) = K + » + b T ' + 2G ) < 82 > 

where we use Tj = (2Tjj_ + Ti||)/3 and G = Ti(Ti\\ — Tjj_)/6 and qq = q^±\\ + q^. The 
anisotropy dissipation term cancels. The first factor of G is from the pressures (diamagnetic 
flow in the fluid model), and the second is from the temperatures (diamagnetic heat fluxes) 
and is neglected in the fluid model. In terms of Pi and G the density equation (Eq. [77)1 
becomes 

df) ■ d 

-± ~ V ■ j^xW + Vpu —JC(<p +Pi + G) (83) 

Subtraction of Eq. ( 1531 from Eq. ( 1201 recovers Eq. (fl9l) above, and then addition of Eq. (fl9l) 
to Eq. (1521) recovers Eq. (1221 above (except for the second factor of G which the fluid model 
doesn't keep), and the correspondence in the nonlinear polarisation terms is thereby proved. 



IV. GYROFLUID TEMPERATURE ANISOTROPY AND FLUID PARALLEL 
VISCOSITY 

We now turn to the less obscure parts of the correspondence between gyrofluid and fluid 
equations. The general pressure tensor arises from 

P = J dWmww /(w) (84) 

where the integration is over velocity space and w is the random kinetic velocity in the 
co-moving frame with fluid velocity u. In the Braginskii fluid equations / is assumed to 
be a Maxwellian f with variable density n and temperature T, and also flow u, where all 
of n, T, u are arbitrarily variable, i.e., including all dynamics as well as the background 
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(reduction to low frequency equations under drift ordering comes later). Then, the correc- 
tions to /o are considered to be of the form fi = $(w)/ , and $ solved for in terms of 
Sonine polynomials, with arbitrarily large collision frequency and small gyroradius, as well 
as small mean free path (i.e., all gradients are assumed to represent small corrections to 
local thermodynamic equilibrium, LTE). Besides the specific conductive heat flux (q/nT), 
the temperature anisotropy is assumed to be small. Hence P is split in terms of an isotropic 
part and a trace-free correction, 

P = pg + IT (85) 

where g is the metric tensor. The diagonal elements of II represent the parallel viscosity. 
These pressure contributions may be written, separately from any non-diagonal II contribu- 
tions, as 

P = P± g + (P\\ - P±)hb (86) 
where b = ~B/ B is the magnetic unit vector, hence 

Ii = Aphh--Apg (87) 
3 

with isotropic pressure and deviation given by 

P 

We also have 

P\\ 

to assist the evaluation of gyrofluid combinations. 

It is important to note that with the inhomogeneous magnetic field V • (bb) contributes 
to the general divergence V • P. Then, in the reduction to low frequency we obtain 

b • (V • P) = V||P|| - ApVy log B (90) 

for the parallel pressure force, which is covered by the parallel gradient and magnetic pump- 
ing terms in the gyrofluid moment equation for the parallel velocity. In the perpendicular 
drifts we make the same approximations as in the gyrofluid model itself: 

V-^BxVp^-2Vlo gj B-^-BxVp b-Vb^Vlog5 (91) 



2p ± + p|| 



Ap = p|| - p ± 



(88) 



p+ -Ap 



P±=P~ g &p 



(89) 
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and hence also 



We can then find 



V--^Bx(b- Vb) -> -2Vlo gj ■ -^-BxVlogE = (92) 
B z B 2 



V • JjBx(V • P) -> -Vlog£ • ^BxV(p|| + Pl _) (93) 

which is covered by the same curvature terms in the gyrofluid moment equation for the 
density. The gyrofluid model then also includes FLR effects, wherein Vy0 becomes Vii^g 
and under derivatives and in the magnetic pumping terms TjT_i_ becomes TiTi± + Qq- 

Hence, the pressure/temperature anisotropy enters the continuity equations for ions the 
same way as in the fluid model, although via a different route: grad-B and curvature drifts for 
gyrocenters rather than the pressure tensor for particles. It remains to obtain the anisotropy 
itself. The gyrofluid equations for Tjj_ and Tjii for the ion species are given by 

lit + ^ [0G ' Tl||] + BW W Ull+ B qiU ~ + g4±l|)V|1 bg B 

r ( <t>G + TiUi + &TiTi\\\ Vi r . 

= [ 2 j ~ 3^' ~ ' ~ G * ^ ' 

FfT n 

+ [<h, Ti±] + [Sla, (ni + 2T U )] + BV,^ + ( M|| + 9j±|| ) V| logB 

= K f *a + Tg; + 4n7U + 4H G \ + ^ _ _ 

\ 2 / 37Tj 

where the Vn log-B terms give the magnetic pumping of anisotropy and the terms its 
collisional dissipation, and 7Tj is a numerical constant which we eventually adjust to obtain 
correspondence. If we add these equations the total temperature equation results, and the 
magnetic pumping and dissipation cancel, leaving 

~± + (FLR) + Wy^-j^ = JC (V + nn t + 7 - n T t + ^At) (96) 

where "FLR" denotes the FLR corrections which eventually become the polarisation terms 
as established above and Tj = (2Tjj_ + Tj||)/3 and qi\\ = qi\\\\ +qi±\\ as before. The anisotropy 
is AT = Tj|| — Tjj_. If we instead subtract Eq. (1951) from twice Eq. f)94p forming an equation 
for AT, we find 

<9AT . . ,. 2q , ) -|iii — 0? 1 11 

— + [<P G , AT] - [Q G , (th + 2T i± )) + BV|| ML_^m 

- 3(«„ + ft± „)V|| io gB - k ^Lz±±Ik _ 1^ _ ^ no 
= \k (0 + nm + nTi) - 2£V,3 - -t,at (97) 
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having arranged terms such that the collisional and velocity divergence terms are on the 
right side and the nonlinearities, magnetic pumping, FLR, anisotropy in curvature terms, 
and heat flux effects are on the left side. The Braginskii assumptions are essentially that the 
left side terms are small, even though that is obviously not the case in the curvature terms, 
as pointed out before }l3j . 

If the Braginskii assumptions are taken, then we find 

(98) 



r . A T = ^ 

2v; 



K (0 + r ini + nTi) - 4BV,3 



which is equivalent to 

7T,: 



G 



(Vlog J B i -u ± )-4(V-M||b) (99) 



This is the form given by Ref. 24J ]. and correspondence is thereby proved. It is important to 
note, however, that this regime is never reached even in deep edge turbulence. Tokamak edge 
turbulence typically has about two orders of magnitude slower than nonlinear advection. 
Worse than this, the nonlinear advection is the largest effect in Eq. (l97l) . larger than any 
of c s /R or c s /qR or z/j, even for zonal flows. Hence the dissipation of ion flows cannot be 
properly modelled by K{G) in the vorticity equation or VmG in the parallel velocity equation, 
with either z/j or c s /qR or c s /R as the controlling frequency. These forms will overestimate 
the dissipation effects and will enforce a particular phase shift between viscosity and the 
variables determining the ion flow (<ft, rij, 7$, and u\\). Use of the collisional form with 
Vi will very strongly overestimate ion flow damping in tokamak edge regimes, ultimately 
corrupting any investigation of bifurcation dynamics. Nevertheless, we set 7Tj = 0.96 for 
ions and 7r e = 0.73 for electrons to obtain correspondence to the Braginskii regime in the 
gyrofluid equations. 



V. GYROFLUID HEAT FLUXES AND THE COLLISIONAL FLUID LIMIT 

The simplest correspondence is in the heat flux equations. In the gyrofluid model the 
parallel heat fluxes (parallel and perpendicular energy components) are dynamical variables 
with their own equations. However, the Braginskii limit assumes that all of the time scales 
involved in advection, divergences, dissipation, etc., of the heat fluxes are slow, with the 
exception of collisional dissipation. The dissipation balances the forcing represented by the 
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temperature gradient. When the parallel heat flux components (Eqs. 103,104 of Ref. 23) are 
added to form the equation for qq = qqu + the anisotropy dissipation effects cancel, 
leaving 

^^If + (FLR) + (LD) ~ ^ r ^ 2u \\ + 4 ^im + 3%± n) 

= — nV|| (rjTq + T i± ^j - ^Hinuiqq (100) 

where Ki is the thermal conduction coefficient, and "FLR" denotes the FLR corrections 
(including the appearance of VLq under Vy) and "LD" the Landau damping dissipation. All 
of the terms on the left side scale with advection or are slower. The terms on the right 
side are the ones left after the Braginskii ordering is taken — assuming that z/j overpowers 
advection. Neglecting the temperature anisotropy the factors of 5/2 cancel and we have 

qq = -— V||T< (101) 

which is the Braginskii formula. We set Ki = 3.9 to set the quantitative correspondence. 

For electrons we additionally have the mixing of the moments under the collisional dis- 
sipation. The re-expression of the thermal force as such was given in Ref. [5l|, whose two 
salient equations in the fluid model are 



„ Mi, dJ\\ „ . m 

Pe^TT + ^e-TT = V||(n e + T e - 0) - [l e V e 

at at 



r)J\\ + — (g e || + a e J\\) 

K P , 



(102) 



da ii 5 5/2 

/i e -^ + (LD) = --V,|T e - J-^q^ + Ml) (103) 
at 2 K e 

where the coefficients 77, K e , a e are for resistivity, thermal conduction, and the thermal force, 
respectively. If the nonlinear advection and Landau damping are assumed small in Eq. ( 11031) 
then the Braginskii formula 

g e |l+« e J|| = -^V|,T e (104) 

with K e = 3.2 and a e = 0.71 for pure hydrogen, is recovered. Then, insertion of this into 
Eq. (jIQ21) gives 

dAii dJ\\ 

Pz-Qf + = V ll( n e + Te - <P) + «eV||T e - T][l e V e .J\\ (105) 

which recovers the Braginskii Ohm's law if rj = 0.51 is chosen. 
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This dissipation model was built into the electron gyrofluid moment equations in order 



to obtain this correspondence, in both Refs. 
dq, 



471 | . The electron heat flux equations are 



Of 



+ fi e a Le q e \\\\ +/i e [0 e ,g e ||||J = --V||i e || - /U e /C I — — — I 



(5/2) 



VeVe (?e|||| + 0.6ct e J||) + 1.28f e (g e |||| - 1.5<? e 



(106) 



<y ' ± " + Ai e a Le g e± || + /i e [0 e ,g e _L||] +yu e [fi e , (w|| +2g e _L||)] 



9* 



-Vn(T e± -n e )-/i e x: 

(5/2) 



uii + 6q e 



(r e± -T e \\ -n e ) Vyiogs 



ft, r .. 



/" e ^e <? e ±ii + 0.4a e Jii - 1.28v e (g e ||i| - l.5q 



(107) 



from Eqs. (103,104) with additions in Eqs. (114,118,119) of Ref. [23]. Electron FLR correc- 
tions are kept, with e and Q e the corresponding potentials. Adding these to form the total, 
neglecting FLR effects, magnetic pumping, and curvature terms, we find 



+ < LD > 



(5/2) 



{J>eVe(q e \\ + OC e J\\) 



(108) 



noting the anisotropy dissipation terms cancel. Now assuming u e overcomes nonlinear ad- 
vection or Landau damping, and neglecting AT, we find 



q e \\ + a e J\\ 



"V||T e 



(109) 



which is the same as Eq. (j!04p above, i.e., the Braginskii formula. 
The electron gyrofluid parallel velocity equation is 



Pe-fo ~ Ve-r^ - He[4>e, V\\] ~ /i e [^ e , q e 



-Vii 



-n e - T e| |) - ix JC 



Av\\ + 2g e i||| + q e 



(p e - T e± + r e ,|) V,| logB-fj, t 



T]J\\ + 



a,. 



+ q e ±\\ + at e J\\ 



;no) 



from Eq. (100) with additions in Eq. (115) of Ref. 23|. Neglecting FLR effects, magnetic 
pumping, curvature terms, adding q e ^ + q e ±\\ = q e \\, replacing T e \\ by T e (neglecting AT) 
and setting uy = —J» (effectively neglecting finite /i e corrections), we find 



„ dA\\ dJ\\ „ , 



a 



VJ\\ H (q e \\ + at e J\\) 



(111) 
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which is the same as Eq. (11021) above. Then going to the Braginskii limit by inserting q e \\ 
from Eq. (11041) . we find Eq. (11051) . which is the Braginskii version. 

It has been pointed out that curvature terms should appear in the fluid model's equations 
for Ju and u\\ [13J. This is indeed the case, but in that event the involved terms are the same 
as the ones in the gyrofluid model. Correspondence in the heat fluxes and the Ohm's law is 
thereby proved. 

VI. CONCLUSIONS 

This work has considered the low frequency Braginskii fluid drift equations on the one 
hand and the electromagnetic, transcollisional gyrofluid equations on the other. In both cases 
the model is comprehensive enough to treat temperature dynamics of both species and has a 
free energy functional which is conserved in the absence of dissipation and external drive. By 
considering the entirety of the Braginskii limit — collision frequency larger than advection 
or transit frequency, specific heat flux smaller than fluid velocity, and small gyroradius - 
the sets of equations and each one's conserved free energy have been shown to be one and 
the same. Whether or not the Braginskii limit is ever reached is a separate question. But 
the low frequency Braginskii fluid drift equations have been shown by correspondence to be 
a fully contained subset of the electromagnetic, transcollisional gyrofluid equations. 

One substantial advantage of the gyrofluid model is the fact that in all the terms involving 
derivatives, only scalar quantities are involved. Instead of vector or tensor components, it 
has gyroreduced potentials or charge densities which involve Hermitian operators (enabling 
the free energy conservation). This leads to numerical schemes which are easier to formulate. 
All the nonlinear terms have the Poisson bracket structure, for which the Arakawa spatial 
discretisation scheme is uniquely suited [52]. Coupled with a timestep that i s hig hly accurate, 



requires only one evaluation of the terms per step, and is stable for waves 53[, we have the 



best scheme found so far for this type of microturbulence in magnetised plasmas 54j. Its use 
for the gyrofluid equations is detailed in Ref. [23j. With these advantages together with the 



correspondence to the Braginskii fluid drift equations, some of the mystery surrounding the 
efficacy or validity of the gyrofluid model for tokamak edge turbulence should be alleviated. 
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